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GRAPHICAL  METHOD  OF  OBTAINING  THEORETICAL 
LIFT  DISTRIBUTIONS  ON  THIN  WINGS 
AT  SUPERSONIC  SPEEDS 
By  Clarence  B.  Cohen  and  John  C.  Eward 


SUMMARY 

A  graphical  method  is  presented  for  calculating  the  linearized 
lift  distribution  on  thin  vings  at  supersonic  speeds.  The  technique 
may  be  applied  to  all  wing  regions  except  those  influenced  by  inter¬ 
acting  flow  fields  off  the  wing-plan  boundaries.  The  lifting- 
pressure  coefficients  are  obtained  as  the  sum  of  a  graphical  line 
integration  and  several  terms  that  are  functions  of  only  the  plan 
form. 


Pressure  coefficients  obtained  by  the  graphical  method  for  a 
swept  trapezoidal  wing  and  for  a  wing  with  straight-swept  leading 
edge  and  parabolic  wing-tip  are  compared  with  those  derived  by 
closed-form  integration.  The  lift  distributions  for  two  plan  forms 
previously  unsolved  are  also  included. 


INTRODUCTION 

The  evaluation  of  the  theoretical  performance  of  arbitrarily 
shaped  wings  at  supersonic  speeds  has  been  retarded  by  the  complex¬ 
ities  involved  in  analysis  when  the  wing  boundaries  are  other  than 
certain  prescribed  curves,  usually  straight  lines.  Solutions  for 
the  velocity  potential  on  an  arbitrary  wing  surface  are  presented 
in  references  1  and  2  except  for  those  regions  influenced  by  inter¬ 
acting  flow  fields  off  the  wing-plan  boundaries;  the  integrals 
involved,  however,  have  been  tabulated  only  for  special  plan  forms. 
Although  the  integrals  may  be  so  reduced  that  numerical  evaluation 
of  the  velocity  potential  is  feasible,  the  procedure  for  obtaining 
the  lift  distribution  requires  an  additional  numerical  partial  dif¬ 
ferentiation  that  is  cumbersome  and  of  doubtful  accuracy. 

The  lift  distributions  for  a  family  of  thin  wings  may  be  obtained 
by  the  methods  of  reference  2  as  the  sum  of  an  algebraic  function 
and  a  line  integral.  The  line  integral  may  be  evaluated  by  numerical 
or  graphical  methods. 
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A  graphical  method  formulated  at  the  NACA  Cleveland  laboratory 
for  evaluating  the  lift  distribution  for  wings  of  arbitrary  plan 
form  within  the  limitations  of  reference  2  is  presented.  Also 
included  for  a  Mach  number  of  *{2  are  curves  showing  the  relations 
between  the  algebraic  functions  and  the  geometry  of  the  wing.  The 
application  of  the  method  to  the  various  types  of  region  encountered 
in  analyses  of  general  wings  is  discussed.  As  illustrations,  the 
pressure  coefficients  obtained  by  this  graphical  method  are  compared 
with  those  derived  by  closed-form  integration  for  a  swept  trapezoidal 
wing  and  for  a  wing  with  a  straight  leading  edge  and  parabolic  wing- 
tip  boundaries.  Lift  distributions  of  two  other  plan  forms  are  also 
included  to  illustrate  the  application  of  the  method. 


SYMBOLS 

The  following  symbols  and  subscripts  are  used  in  this  report: 

a2>  a3>  •  •  •>  an  coefficients  in  series  expansion  defining  curva¬ 
ture  of  leading  edge 

a',  b',  c',  d'  quantities  used  in  determining  factor  AF2 

Ci  .  .  .  Cn  coefficients  in  series  representation  of  AF2 

Cp  pressure  coefficient 

AGp  contribution  to  Cp  of  portion  of  leading  edge 

included  in  segment  of  width  e 

factor  representing  integral-function  contribu¬ 
tion  to  Cp 

Fg  factor  representing  contribution  of  integration 

over  width  e  to  Cp 

AF 2  error  in  Fg  due  to  neglecting  curvature 

H  quantity  representing  numerical  integration 

along  leading  edge, .  ArjZh 

h  parameter  defining  equilateral  hyperbola 

kq  constant  determining  sweep  of  straight  leading 

edge 


972 


NACA  TN  Wo.  1676 


3 


M 

m 

Q(c) 

H 

r 

U 

u,  v 

u2(v)  or  v2(u) 
v^(u)  or  u^(v) 

x,  y 

r 

a 

0 

€ 

6l>  e2 

n 

cp 

Subscripts: 

0 

L 


free-stream  Jfech  number 

slope  of  wing  edge,  ^ 

series  defining  curvature  of  leading  edge 

vw  -  V-jCug) 

ratio  of  distances,  - * 

'  %  “  «2 

local  radius  of  curvature 

free-stream  velocity  (parallel  to  x-axis) 

oblique  coordinates  whose  axes  lie  parallel  to 
Mach  lines 

equation  of  tip  and  trailing  edge  of  wing 

equation  of  supersonic  leading  edge  of  wing 
(inboard  of  foremost  tangent  Mach  wave) 

Cartesian  coordinates 

transformed  y- coordinate 

angle  of  attack,  radians 

Vm2  -  1 

width  of  strip  chosen  in  evaluating  F2 

local  angles  between  wing  edge  and  free-stream 
flow  direction 

Cartesian  coordinates  of  point  sources 
perturbation- velocity  potential 

intersection  of  grid  edge  with  wing-plan  boundary 
leading  edge  of  wing 
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l 

left 

r 

right 

w 

wing 

ANALYSIS 

Aa  a  result  of  the  linearized  theory,  the  effect  of  angle  of 
attack  on  the  lift  distribution  of  thin  wings  at  supersonic  speeds 
depends  upon  the  plan-form  boundaries  and  is  independent  of  the 
wing-section  slopes.  Determination  of  the  pressure  distribution  of 
a  thin  flat  plate  therefore  solves  the  problem  of  the  lift  distri¬ 
bution  of  the  finite- thickness  wing.  Eeference  2  shows  that  the 
local  pressure  coefficients  of  several  classes  of  wings  with  arbi¬ 
trary  plan  boundaries  may  be  evaluated  by  line  integrals  of  explicit 
functions . 

A  graphical  method  of  evaluating  the  lift-distribution  functions 
presented  in  reference  2  is  developed  herein.  In  order  to  illus¬ 
trate  the  method  developed,  a  simple  wing  including  regions  influ¬ 
enced  by  a  supersonic  leading  edge,  a  subsonic  leading  edge,  and  a 
subsonic  trailing  edge  (fig.  l)  is  discussed  in  detail.  The  essen¬ 
tial  elements  in  the  calculations  for  each  region  are  included  in 
the  analysis  for  the  region  B  influenced  by  the  subsonic  leading 
edge. 


The  pressure  coefficient  for  region  B  as  derived  in  refer¬ 
ence  2  is 


r  _  Ma 
°P  ~  ”  afi 


L 


dji 


ab  V(*v  -  u)  (vw  -  v) 


(1) 


where  ug  is  evaluated  at  v  =  vw.  Except  for  the  symbol  tj, 
equation  (l)  is  expressed  in  a  set  of  oblique  coordinates  (fig*  2) 
whose  axes  are  parallel  to  the  Mach  lines*  The  transformation 
equations  relating  the  oblique  and  Cartesian  coordinates  are 
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The  significance  of  the  terms  of  equation  (l)  is  clarified  in 
figure  3.  The  line  integral  is  to  be  evaluated  along  the  portion 
of  the  supersonic  leading  edge  from  a  to  b.  The  quantity 
[vw  -  vq ( ug ) ]  can  be  interpreted  as  distance  be  and  the  quantity 
(uy.  -  Ug)  can  be  interpreted  as  the  distance  cP.  The  ratio  of 
these  distances  is  called  E.  For  a  given  Mach  number,  the  quan- 

f  du2\ 

tity  f  1  -  rj— -  )  depends  upon  only  the  slope  of  the  wing  boundary 
at  the  point  (u2,vw). 


By  use  of  equations  (2),  the  second  term  of  equation  (l)  may 
be  expressed  in  terms  of  E  and  the  angle  0q  that  the  wing  edge 
makes  with  the  free-stream  flow  direction  (parallel  to  the  x-axis) 

as 


4p_tan_ei_f- 

=  M  (p  tan  0q  +  l)11 


(3) 


where  0q  is  positive  counterclockwise.  This  expression,  eval¬ 
uated  for  M  =  v2,  appears  in  figure  4. 

The  quantities  (uy.  -  u)  and  (vv  -  v)  in  the  line  integral 
of  equation  (l)  are  the  coordinate  distances  (fig.  3)  of  each 
element  dq  to  the  point  (uw.,vw).  The  curves  of 


V(« w  "  u)  (vw  -  v) 


=  constant  =  h 
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in  the  u,v  coordinate  system  are  a  family  of  equilateral  hyper¬ 
bolas  asymptotic  to  the  lines  u  =  u^,  v  =  vw.  A  family  of  these 
hyperbolas  is  shown  in  figure  5. 

If  a  "grid"  of  this  type  is  placed  with  its  origin  at  (uH->vw) 
and  so  alined  that  its  axes  are  parallel  to  the  u  and  v  axes, 

the  value  of  — - — -  for  any  element  dq  (or  dy)  is 

VW  -  u)  (vw  -  y) 

the  constant  h  for  the  hyperbola  passing  through  that  element. 

The  sum  of  the  values  of  b  taken  near  the  center  of  each  element 
for  equal  increments  drj  along  the  wing  boundary  then  evaluates 
the  line  integral  of  equation  (l).  The  superposition  of  a  wing- 
plan  boundary  on  a  hyperbolic  grid  for  M  =  V2  is  shown  in  fig¬ 
ure  6.  A  constant  value  of  Atj  is  marked  on  the  grid  for  conven¬ 
ience  in  summation.  A  part  of  the  grid  employed  for  calculations 
at  M  =  V2  appears  in  figure  7 . 

As  shown  in  figure  5,  the  value  of  h  increases  to  infinity 
as  the  edge  of  the  grid  is  approached.  In  order  to  perform  a 
numerical  integration,  the  infinities  must  be  excluded.  Therefore 
the  numerical  summation  should  not  start  with  element  1  (fig.  6) 
but  should  skip  a  distance  large  enough  to  enable  the  element 
value  to  be  easily  determined.  The  part  of  the  integral  omitted 
must  be  compensated  by  an  analytical  expression  that  represents 
the  contribution  of  a  section  of  the  leading  edge  aa*  included 
in  a  strip  of  width  e  (fig.  6). 

The  contribution  to  the  line  integral  of  the  strip  of  width  e 
may  be  evaluated  by  assuming  the  wing  boundary  to  be  a  straight 
line  near  u  =  u^.  From  the  derivation  presented  in  appendix  A, 
this  contribution  is 


Fo  = 


4  pm 


M 


V(m0  -  l)(mp  +  1) 


tan' 


-1 


2|3m(vw  -  t0) 

M  €  (1  +  pm) 


1 


dy 

where  m  =  g,  the  tangent  of  the  angle  ©3  (positive  counter¬ 
clockwise)  that  the  leading  edge  makes  with  the  flow  direction  at 
the  intersection  with  the  edge  of  the  grid  (point  a,  fig.  6),  and 
(vw  -  Vq)  is  the  distance  aP  from  that  point  to  the  grid  origin 

(fig.  6).  The  value  of  F2  for  various  values  of  ©3*  (vw  *  v0). 
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and  e  Is  shown  in  figure  8  for  M  =  V2.  The  dashed  line  in  fig¬ 
ures  8(h)  and  8(c)  represents  the  limit  for  which  the  arc- tangent 
term  in  equation  (4)  is  real. 

For  the  subsonic  leading  edge  (region  B,  fig.  1),  equations  (l), 
(3),  and  (4)  may  he  combined  and  the  expression  for  the  pressure 
coefficient  becomes 


Cp  =  -  ^  (H  +  FL  +  F2)  (5) 

where  the  numerical  summation  along  the  leading  edge  is  repre¬ 
sented  by 

H  =  At)  (hp  +  h2  +  h3  +  .  .  .)  =  ArjSh 

If  nonunifora  increments  of  t)  are  desired,  a  weighting  of  each 
value  of  h  is  required. 

The  analysis  of  the  region  of  a  wing  tip  influenced  by  a  sub¬ 
sonic  trailing  edge  (region  C,  fig.  l)  may  be  appreciably  shortened 
if  the  Kutta-Joukowski  condition  is  imposed.  In  reference  2,  it 
is  shown  that  the  solution  which  satisfies  the  Kutta-Joukowski  con¬ 
dition  does  not  contain  the  tera  Fp.  Therefore  Fp  is  taken  as 
zero  in  this  region  and  the  remaining  evaluation  is  unchanged. 

When  evaluating  Cp  in  a  region  where  v<  0  (region  A, 
fig.  1,  or  fig.  9),  the  value  Fp  does  not  exist  because  no  sub¬ 
sonic  leading  edge  is  included  in  the  f orward  Mach  cone  from 
point  P.  The  integration  must  be  conducted  over  the  entire  wing 
leading  edge  between  the  limits  of  the  grid  a  and  b.  Because 
h  becomes  infinite  at  the  right  as  well  as  at  the  left  limit  of 
the  grid,  a  part  of  the  integration  at  the  right  must  be  replaced 
by  an  analytical  expression  similar  to  equation  (4).  Equation  (4) 
can  be  shown  to  hold  at  this  limit  if  (uy-  -  uq)  is  substituted 
for  (vw  -  vq)  and  (180°  -  8p)  is  substituted  for  02.  If  the 
value  of  this  expression  is  designated  F2  r  the  pressure  coef¬ 
ficient  in  the  region  where  v<0  becomes 

Cp  =  -  ^  (H  +  F2  +  F2^r)  (5a) 

In  general,  a  complete  wing  In  supersonic  flow  may  have  regions 
under  the  influence  of  both  wing  tips.  A  summary  of  the  method  of 
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obtaining  lift  distributions  in  the  various  types  of  flow  field 
commonly  encountered  is  given  in  appendix  B.  Also  included  is  a 
numerical  evaluation  of  the  pressure  coefficient  at  point  P  of 
figure  10. 

The  effect  of  Mach  number  on  lift  distribution  for  a  given 
plan  form  may  be  determined  by  two  alternate  procedures.  One 
method  is  to  construct  an  integration  grid  and  curves  of  Fp  and 
Fg  for  the  Mach  number  desired.  If  a  large  number  of  analyses  at 
a  given  Mach  number  are  to  be  made,  this  method  may  be  preferable. 
If  only  a  very  few  analyses  at  a  given  Mach  number  are  to  be  made, 
however,  the  grid  and  curves  of  Fp  and  Fg  may  be  used  as  con¬ 
structed  for  M  =  V2,  provided  that  compensating  corrections  are 
made.  These  corrections  consist  in:  (l)  so  transforming  the  wing 
boundary  that  y’  =  3y;  and  (2)  analyzing  the  wing  as  at  M  =  V2 
and  dividing  the  resultant  Cp/a  by  the  value  of  3. 

Some  considerations'  that  simplify  the  application  of  the 
method  appear  in  appendix  C.  Construction  of  the  integration  grid, 
choice  of  At]  and  e,  drawing  size,  and  a  special  example  of  the 
subsonic  trailing  edge  are  discussed. 


ACCURACY  OF  METHOD 

In  order  to  determine  the  accuracy  of  the  method,  the  pressure- 
coefficient  at  points  on  two  wings  for  which  analytical  expressions 
were  obtainable  were  graphically  computed  and  compared  with  the 
analytical  values.  The  wings  and  the  points  considered  appear  in 
figure  11.  The  leading  edges  of  both  wings  are  swept  back  30°  and 
are  the  same  except  that  wing  A  has  a  straight  tip  and  wing  B  has 
a  parabolic  tip  passing  through  similar  points  on  the  leading  and 
trailing  edges.  The  points  considered  are  at  the  same  x,y  values 
for  each  wing.  The  results  of  the  computations  for  Cp  are  pre¬ 
sented  in  table  I. 

The  magnitude  of  the  personal  error  incurred  in  application 
of  the  method  was  evaluated  by  trial.  Values  of  Cp/a  were  inde¬ 
pendently  obtained  by  nonprofessional  personnel  using  the  graphical 
method,  and  were  also  analytically  calculated.  For  wing  A,  the  results 
of  one  of  the  computers  contained  an  average  error  of  0.46  percent 
and  a  maximum  error  of  0.71  percent  (table  I).  Another  computer 
averaged  an  error  of  0.53  percent  with  a  maximum  error  of  1.52  per¬ 
cent.  Results  of  wing  B  gave  about  the  same  accuracy  (table  I). 

These  results  were  obtained  from  drawings  with  a  wing  chord  of 
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4g  inches.  The  computers  estimated  a  time  rate  of  approximately 

5  minutes  per  point  when  a  large  number  of  points  were  being 
analyzed. 

Inasmuch  as  the  f irst-order  expression  for  F2  (which  also 
includes  higher-order  terms)  assumed  a  straight  leading  edge,  this 
expression  was  exact  for  the  wings  shown  in  figure  11.  The  accuracy 
of  this  expression  for  a  curved  edge  has  been  investigated  and  is 
shown  in  appendix  A.  For  a  parabolic  leading  edge,  the  error  AF2 
in  Fg  (including  fourth-order  terms)  as  a  function  of  the  ratio 
e/rQ  (where  tq  is  the  radius  of  curvature)  is  presented  In  fig¬ 
ure  12.  If  c/tq  is  small,  the  error  in  Fg  due  to  curvature 
will  be  small  except  when  6%  approaches  the  Mach  angle,  as  indi¬ 
cated  in  figure  12.  The  corresponding  error  in  Cp  will,  of  course, 
be  a  smaller  percentage.  A  similar  trend  probably  holds  for  edges 
of  different  curvature. 


APPLICATIONS 

As  examples  of  the  method,  the  lift  distributions  2Cp/a  of  a 
circular-plan-form  wing  and  a  straight  wing  with  a  circular  tip  are 
shown  in  figures  13  and  14,  respectively,  at  M  =  ^2. 

For  the  circular  wing,  in  the  region  inboard  of  station  D,  the 
pressure  at  a  given  chordwise  station  increases  in  the  outboard 
direction.  For  a  given  spanwise  station,  the  pressure  is  a  maximum 
at  the  leading  edge  and  drops  off  in  the  flow  direction.  The  rate 
of  this  decrease  is  a  minimum  at  the  wing  root. 

Lines  of  constant  pressure  for  the  circular-tip  wing  swept 
back  30°  at  a  Mach  number  of  V2  are  shown  in  figure  14.  The 
Kutta-Joukowski  condition  was  assumed  for  analysis  of  regions 
influenced  by  the  subsonic  trailing  edge.  The  only  region  of  high 
pressures  and  high  pressure  gradients  is  a  small  leading-edge 
region  between  the  tangency  of  the  foremost  Mach  wave  and  the  start 
of  the  subsonic  trailing  edge.  The  pressure  gradient  over  most  of 
the  tip  is  roughly  in  the  spanwise  direction. 


CONCLUDING  DISCUSSION 

A  graphical  method  based  on  linearized  supersonic-flow  theory 
has  been  developed  for  calculating  the  lift  distributions  on  thin 
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wings  at  supersonic  speeds.  The  method  is  applicable  to  all  of  the 
regions  of  arbitrary  wings  except  those  affected  by  interacting 
flow  fields  off  the  wing-plan  boundaries.  The  determination  of 
the  lifting-pressure  coefficient  for  a  given  point  has  been  found 
to  require  5  to  10  minutes  with  a  resultant  average  error  of  less 
than  1  percent  as  determined  by  comparison  with  known  analytical 
solutions. 


Flight  Propulsion  Research  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  May  13,  1948. 
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APPENDIX  A 

CONTRIBUTION  TO  Cp  OF  LEADING-EDGE  SEGMENT 
INCLUDED  IN  STRIP  OF  WIDTH  e 

In  order  to  determine  the  contribution  to  Cp  of  a  segment  of 
width  e,  the  integral  of  equation  (l)  must  be  evaluated.  In  the 
x,y  coordinate  system,  this  integral  is 


The  integral  is  considered  in  the  region  near 

(x  -  |0)2  -32  (y  -  tj0)2  =  0  (A2) 

The  subscript  0  refers  to  the  point  of  intersection  of  the  grid 
edge  and  the  wing-plan  boundary.  In  particular,  the  positive  root 
that  defines  the  left  forward  Mach  line  from  the  point  (x,y)  is 

(x  -  |0)  =  0  (y  -  tj0)  =  |  (vw  -  v0)  (A3) 

Treatment  of  e  as  a  variable  in  the  vicinity  of  the  point 
(t)0>  t  o )  yields 

n  =  ^0  +  e 

|  =  lo  +  |  +  Q(0 

where 

Q(c)  =  age2  +  a3c3  +  a4  e4  +  .  .  . 


and 


12 


NACA  TIT  Ho.  1676 


Substitution  of  equation  (A4)  in  equation  (Al)  gives 


How,  if  Q(  e)  =  0,  the  leading  edge  is  straight  and  the  inte¬ 
gral  may  be  written 


This  function,  which  is  equation  (4)  of  the  text,  is  shown  in  fig 
ure  8  for  various  values  of  and  (vw  -  Vq) / €  for  M  =  v2. 

If  the  leading  edge  is  curved,  Q(e)  4  0  and  equation  (A5) 
must  be  rearranged.  The  quantity  in  the  radical  may  be  expressed 
as 
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where 


a'  =  (o'  +  a*  e) c  = 


LL: 


2p2  (y  -  tj0)  - 


2  (x  -  |0) 


m 


(]  fJL 
j  +  (m2  " 


p2)€ 


>  € 


b.=2U-£0)4-2(|)Q- 
2  (x  -  | 0) (mp  -  1) 


c  *  = 


,  _  (1  -  nr  p^) 


m 

2  n2 


d'  = 


(A8) 


m1 


Then  equation  (AS)  may  be  written 


Expanded  in  a  power  series, 


Equation  (A5)  then  becomes 


The  first  term  is  identical  with  the  integral  of  equation  (A6). 

The  contribution  of  the  curvature  to  F2  may  therefore  be  written 
as 
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The  expressions  for 


V 

a' 


expanded  in  a  power  series  to  yield 


in  equation  (AS)  may  be 


_2_ 

o' 


[(x  -  |0)a2]e  +  (x  -  l0)a3  “  ~  "  “7  <x  "  ^  o)a2  J  * 


+  [(*-£  0)&4  "IT  a22  -  (fO  (x  -  &0)a3  +  (fj)  ^ 

+  (£r)2  (x-  + 

_i _ i _ i _ 

^  ^  f  (1  4  41  C) 


_1  _1_ 

Vc  Vc^ 


where 


41  =  -d  ±  mg)  (A10) 

c'  2m  (x  -  |Q) 

Upon  substitution  of  equation  (A10)  in  equation  (A9)  and  multi¬ 
plication  and  collection  of  terms  according  to  powers  of  €,  there 
results 


?.)A 
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* 


3  rt 


^2  =  if  (-f)2 J0  (3Cle2  +  5C2e2  +  7C3e 


+  * 


d€ 


—  /  KS  SJ  * 

—  {.  I  p  ,p(r2.p^2, 

"m  Uv  V1  °2  °3 

where 

°1  =  k  (x  “  lo^a2 


°2  =  i 


C3  =  1^ 


(36m  -  l)a2 


3a?2  id 

Ct 

a*  + - 

(x  -  |0j  ► 

4 (pm  -  1)_ 

(A9a) 


4m 


3(|3m  +  l)(50m  -  3)a2  (7pm  -  l)a2  (3pm  -  l)a3 

—  +  — : -  +  - 


32m2(x  -  |0) 


16(pm  -  1) 


4m 


a4  + 


^a2a3la 


3  2 


5a  o  m 


2(pm  -  1)  8(pm  _  X)2  J 


(x  -  lo) 


(All) 


C4  = 


The  series  of  equation  (A9a)  converges  as  long  as  b'/a*  and. 
d'e/c'  (defined  by  equations  (AS)  and  (A10))  are  less  than  unity, 
For  rapid  convergence,  these  ratios  should  both  be  small. 

As  a  particular  example  (taken  at  M  =  V2),  the  leading-edge 
curve  may  be  represented  by  the  parabola 


1=  £q  +  £  +  a2e2 
m 

The  coefficients  of  equation  (All)  evaluated  in  terms  of  the  oblique 
coordinate  distances  then  become 
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f 


J 2 

C1  “  T  *VV  “  T0^a2 


5  V2  a23m2(yw  -  vQ) 
16 (m  -  l)2 


When  the  values  of  the  coefficients  are  substituted  In  equation  (A9a ) , 
AF2  results. 

For  a  parabola  with  a  given  initial  slope  m  at  the  point 
(!qj  t|q) ,  the  quantity  a2  controls  the  radius  of  curvature  Vq 
at  that  point.  The  ratio  AF2/f2  represents  the  relative  error 
resulting  from  the  assumption  of  a  straight  leading  edge  for  a 
width  0  when  the  actual  edge  is  curved.  This  value  from  equa¬ 
tions  (A7)  and  (A9a)  has  been  plotted  against  the  ratio  e/rQ  in 
figure  12. 


'2  =  5 


(3m  -  l)a2  3  a^fvy.  -  vQ) 
-  + - 


4m 


8(m  -  1) 


C3  “  7 


3 +  l)(£m  -  3)ag  (7m  -  ljag2 

-  +  -  H 

32m2  (vw  -  vQ)  16 (m  -  1) 
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APPENDIX  B 


■TREATMENT  OF  COMPLETE  WINGS 

A  tabular  summation  of  the  various  flow  fields  encountered  in 
wing  problems  is  presented  in  table  II  and  the  graphical  solution 
in  each  region  is  indicated.  In  the  reference  sketches  for 
M  =  42,  P  is  the  point  at  which  the  pressure  coefficient  is  to  be 
evaluated.  The  numerical  line  integration  along  the  supersonic 
leading  edge  is  performed  from  point  a  to  point  b.  The  direc¬ 
tion  of  integration  of  the  line  integrals  is  indicated  by  arrow¬ 
heads  on  heavy  lines.  For  a  subsonic  trailing  edge  (sketches  C,  F, 

G,  H),  the  ICutta-Joukowski  condition  was  imposed  and  made  F]_  =  0 
for  that  edge.  The  factor  F2  vanishes  when  the  limit  of  the 
leading-edge  line  integral  is  other  than  the  integration-grid  edge 
(sketches  D,  E,  F,  G,  H) ,  because  no  infinities  are  then  encountered. 
When  the  forward-reflected  Mach  lines  from  the  point  P  cross  on 
the  wing  surface  (sketches  E,  F),  the  integration  H  is  considered 
negative. 

The  occurrence  of  these  flow  fields  on  a  schematic  wing  is 
illustrated  in  figure  15.  Each  flow  region  is  identified  by  the 
corresponding  letter  from  table  II.  In  this  example,  flow 
region  H  does  not  occur. 

As  an  example  of  the  method,  nurnerJcal  calculations  of  the 
pressure  coefficient  at  point  P  for  the  wing  boundary  shown  in 
figure  10  are 

01  =  7.5° 
be  =  0.0419 
cd  =  0.084 


H  =  ==  =  0.4988 
cd 
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92  =  125.5° 

(TV  -  v0)  =  0.0772 

e  =  0.01 

H  =  0.01  [(0.4) (31)  +  25.0  -t-  21.0  +  18.85  +  17.55 
+  16.75  +  16.50  +  16.6]  =  1.446 

_f£  =  +^(H+F1+  F?)=  1.043 

art  A 


p 
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APPENDIX  C 


SIMPLIFYING  CONSIDERATIONS  IN  APPLICATION  OF  METHOD 

A  simple  method  of  obtaining  coordinates  of  hyperbolas  to  con¬ 
struct  an  integration  grid  employs  log- log  graph  paper.  If  a 
straight  line  of  slope  equal  to  -1  is  drawn  on  this  paper,  the 
coordinates  of  this  line,  when  plotted  in  Cartesian  coordinates, 
describe  the  desired  equi lateral  hyperbola. 

The  addition  to  the  grid  of  lines  of  constant  v  and  u 
(parallel  to  the  grid  edges)  with  accompanying  scales  at  the  grid 
edges  increases  the  rapidity  of  evaluation,  because  all  signifi¬ 
cant  distances  may  then  be  read  directly.  A  part  of  the  grid 
employed  is  given  in  figure  7. 

In  choosing  the  value  of  e ,  the  first  consideration  is  that 
€  be  large  enough  that  the  value  of  h  may  be  easily  determined 
for  the  first  subsequent  element.  It  is  also  seen  from  figure  8 
that  if  e  is  taken  aa  a  simple  decimal  such  as  0.01  or  0.02,  and 
(vw  "  v0)  is  known,  the  ratio  (vw  -  vq)A  be  rapidly  computed 
and  hence  Fg  is  rapidly  determined.  From  figure  6,  point  a 
will  not,  in  general,  fall  at  an  integral  multiple  of  Aq;  then, 
if  €  =  NAq  where  N  is  an  integer,  point  a'  will  not,  in 
general,  fall  on  an  integral  multiple  of  Aq.  Thus  a  section  whose 
width  is  less  than  Aq  appears  before  the  first  complete  leading- 
edge  element  and  remains  to  be  included  in  the  numerical  integra¬ 
tion.  If  this  small  element  is  evaluated  in  the  normal  manner  and 
weighted  according  to  its  width,  no  discrepancy  will  occur.  This 
weighting  process  is  illustrated  by  the  example  presented  in 
appendix  B. 

The  most  prevalent  personal  error  in  the  application  of  the 
method  was  in  the  measurement  of  angles  and  because 

on  small  drawings  it  is  difficult  to  Judge  accurately  the  tangency 
of  a  curve.  Large-scale  drawings  are  therefore  recommended  for 
measuring  the  angles.  These  drawings  may  also  be  used  for  the 
integration,  although  the  increase  in  accuracy  may  not  Justify  the 
accompanying  inconvenience. 

For  a  straight  supersonic  leading  edge,  equation  (1)  is 
integrable  and  the  expression  for  Cp  becomes 
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du2\  /tw  +  kLu2  +  1  _x  -  u2) 

dTw/  V  %  “  u2  +  an  V  Tw  +  klU2 


(Cl) 


In  the  region  affected  hy  the  subsonic  trailing  edge  (region  abc 
of  'figs.  16  and  17)  the  first  term  in  equation  (Cl)  vanishes  if 
the  Kutta-Joukowski  condition  is  imposed.  The  quantity  in  the 
radical  of  the  arc-tangent  term  is  proportional  to  the  ratio  of 
distances  sP/rs  (fig.  16)  and  hence  any  wing  composed  of  the 
same  leading  edge  and  a  straight  line  through  point  s  would  give 
the  same  value  of  Cp  for  any  point  on  line  sP.  If  lines  of 

constant  pressure  are  desired,  the  arc-tangent  term  in  equa¬ 
tion  (Cl)  must  be  differentiated  and  set  equal  to  zero.  The  result 
involves  a  derivative  du2/dvw,  which  is  determined  by  the  slope 
of  the  edge  at  point  s.  Equivalent  lines  of  constant  pressure 
along  line  s?  may  thus  be  obtained  by  replacing  the  tip  by  a 
straight-line  tip  that  is  tangent  at  point  s.  The  pressure  along 
any  v  =  constant  line  is  therefore  conical  about  a  point  deter¬ 
mined  by  the  intersection  of  the  extension  of  the  leading  edge  and 
the  tangent  to  the  tip  at  n2(vw).  Along  v  =  Vp  and  v  =  Vg 
(fig.  17),  the  constant-pressure  lines  are  conical  about  points  tp 
and  tg.  By  this  relation,  the  constant-pressure  lines  may  be 
quickly  constructed  and  evaluation  is  necessary  only  along  one  line 
of  constant  x  to  determine  the  pressure  field.  The  relation  holds 
only  in  the  region  influenced  by  the  subsonic  trailing  edge  when 
the  Kutta-Joukowski  condition  applies,  for  differentiation  of 
equation  (Cl)  with  the  first  term  included  yields  a  second  deriv¬ 
ative  that  generally  cannot  be  evaluated  from  straight-line 
relations. 
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TABLE  I 

VALUES  OF  Cp/a  OBTAINED  FOR  WINGS  A  AND  B 


Station 

u 

Wing  A 

Wing  B 

Cp/a 

Error 

(percent) 

Cp/a 

Error 

(percent) 

Analytic 

Graphic 

Analytic 

Graphic 

1 

8 

2.4570 

2.468 

0.44 

2.457 

2.447 

0.16 

2 

7 

1 

2.1354 

2.148 

.59 

1.866 

1.875 

.48 

3 

6 

2 

1.9826 

1.975 

.38 

1.674 

1.679 

.30 

4 

5 

3 

1.8577 

1.871 

.71 

1.578 

1.584 

.38 

5 

4 

4 

1.7861 

1.775 

.62 

1.744 

1.733 

.63 

6 

3 

5 

2.3877 

2.400 

.52 

7 

8 

2 

2.0395 

2.033 

.32 

1.775 

1.759 

.90 

8 

7 

3 

1.9295 

1.933 

.18 

1.689 

1.682 

.41 

9 

6 

4 

1.8367 

1.827 

.53 

1.676 

1.670 

.36 

10 

5 

5 

1.7861 

1.797 

.61 

1.862 

1.865 

'  .16 

11 

4 

6 

2.0147 

2.012 

.13 

Average 

■ 

■ 

0.46 

0.42 
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TABLE  II  -  PRESSURE  COEFFICIENT  FOR  EACH  WING  REGION 


Typo  of 
region 

Reference  figure 

Pressure  coefficient 

Special  features 

A 

CP=  §[H+P2^2,^ 

Affected  only  by  supersonic 
leading  edge. 

B 

cp=^[h  +  p1+p2] 

Affected  by  subsonic  and 
supersonic  leading 
edges. 

C 

-  8^—^^ 

_ 

CP  =  MH  +  F2] 

Affected  by  subsonic  trail¬ 
ing  and  supersonic  lead¬ 
ing  edges.  Kutta- 
Joukowskl  condition 
imposed.  F^  =  0. 

D 

CP=  ^[H  +  Pl  +  Fl'll 

Affected  by  supersonic  and 
both  subsonic  leading 
edges,  _F^  j  found  for 
strip  ad  in  same  manner 
as  F-,  for  strip  Be, 
aH 

For  Ft  7  use  R7  =  = 

Ipt  4  dP 

and  measure  0^  j  clock¬ 
wise  at  point  d. 

E 

cp=  ^[-h+f1  +  Fi>1] 

Affected  by  supersonic  and 
both  subsonic  leading 
edges,  F^  j  determined 
as  for  region  D  with 

R  -  I® 

Ri  IF* 

972 


TABLE  II  -  PRESSURE  COEFFICIENT  FOR  EACH  WING  REGION  -  Concluded 


Type  of 
region 

Reference  figure  1 

Pressure  coefficient 

Special  features 

F 

S=  ^[-H  +  Fi,»] 

Affected  by  one  subsonic 
leading  and  one  sub¬ 
sonic  trailing  edge. 
i  determined  as  _ 

previously  with  Rj  = 

Kutta-Jowkowski  con¬ 
dition  imposed  on  sub¬ 
sonic  trailing  edge. 

Fx  =  0. 

G 

°„ .  &  [-] 

Affected  by  both  subsonic 
trailing  edges.  Kutta- 
Joukowski  condition 
Imposed. 

P1  =  P1 ,1  =  0- 

H 

CP  =  w  [H] 

Affected  by  both  subsonic 
trailing  edges.  Kutta- 
Joukowski  condition 
imposed. 

Fi  =  pi,i  -  °* 

I 

Shaded  region  Is  affected 
by  externally  Interact¬ 
ing  flow  fields. 

Solution  not  handled 
by  graphical  method. 
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Figure  1.  -  Illustration  of  wing  regions 
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i 


* 


Figure  3. 


-  Geometric  interpretation  of  terms  in  equation  (1). 
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Figure  7.  -  Part  of  Integration  grid  employed  for  calculations  at  M  mVz. 
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(c)  Range  of  F2  from  0  to  8.8. 

Figure  8.  -  Concluded.  Variation  of  factor  F2  with  (vw-v0)/e  for  various 
values  of  ©g  at  M  :  4i. 
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Figure  10.  -  Determination  of  Cp  for  specific  wing  boundary 

at  M  * 


AVP2 
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Figure  13.  -  Lift  distribution  on  thin  circular-plan-form  wing  at 
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Figure  14 . Lift  distribution  for  M 
2— £  on  subsonic  leading  edge 


•V»  on  circular- tip  wing  swept  back  at  30°. 
=  oo,  on  subsonic  trailing  edge  s  0. 
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Figure  15.  -  Wing  illustrating  regions  considered  in  table  II 


